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Suppressing the lower Multipoles in the CMB Anisotropies 
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The Cosmic Microwave Background (CMB) anisotropy power on the largest angular scales ob- 
served both by WMAP and COBE DMR appears to be lower than the one predicted by the standard 
model of cosmology with almost scale free primordial perturbations arising from a period of infla- 
tion One can either interpret this as a manifestation of cosmic variance or as a physical 
effect that requires an explanation. We discuss various mechanisms that could be responsible for the 
suppression of such low £ multipoles. Features in the late time evolution of metric fluctuations may 
do this via the integral Sachs- Wolfe effect. Another possibility is a suppression of power at large 
scales in the primordial spectrum induced by a fast rolling stage in the evolution of the inflaton 
field at the beginning of the last 65 e-folds of inflation. We illustrate this effect in a simple model of 
inflation and fit the resulting CMB spectrum to the observed temperature-temperature (TT) power 
spectrum. We find that the WMAP observations suggest a cutoff at k c = 4.9j^J x 10 _4 Mpc _1 
at 68% confidence, while only an upper limit of k c < 7.4 x 10 _4 Mpc _1 at 95%. Thus, although 
it improves the fit of the data, the presence of a cutoff in power spectrum is only required at a 
level close to 2a. This is obtained with a prior which corresponds to equal distribution wrt fc c . We 
discuss how other choices (such as an equal distribution wrt In k c which is natural in the context of 
inflation) can affect the statistical interpretation. 
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I. INTRODUCTION 

The most appealing cosmological scenario emerged in 
the mid 1980s as the flat cold dark matter model with 
a nearly scale free initial power spectrum, characterized 
by a gravitational potential £; 3 |<I>fc| 2 = P$(/c) = A s fc™ a_1 , 
with n s « 1. Besides the usual ingredients like baryon 
content £lg and a present expansion rate Hq, the model 
has only a single free parameter, namely the amplitude 
of initial scalar perturbations A s . However, the wealth of 
observational results over the past decade has forced the 
inclusion of an extra parameter in the form of an effective 
cosmological constant 

Inflation plays a crucial role in explaining the flatness 
and homogeneity required by the observations, as well 
as providing a theory of theprimordial spectrum of cos- 
mological perturbations HIlHSSIESEJ 
The simplest single field models of inflation 
predict a flat universe with primordial perturbations 
which slightly deviate from exact scale invariance, n s « 
0.95 ± 0.03 or so while more elaborate models with 
more parameters can go beyond this interval. For in- 
stance in the simplest models of hybrid inflation with 
two fields one can have n s > 1 16]. A long period of 
inflation (with more than 65 e-folds of expansion) gener- 
ates a huge patch which is much larger than the presently 
observable universe. 




CMB experiments manifest a strikin g a greement with 
this standard cosmological model However, an 

interesting result of the WMAP data, which confirms ear- 
lier COBE DMR observations, is a lower amount of power 
on the largest scales when compared to that predicted by 
the standard ACDM models Q,HI3. Accurate Monte 
Carlo simulations of the WMAP observations indicate 
that only 0.7% of the realizations of the models studied 
in |3j have less power than the observed quadrupole. This 
number, which approximates the probability of observing 
a smaller quadrupole than the one observed by WMAP 
in a realization of standard ACDM model, should not be 
compared to the posterior probability for the theoretical 
quadrupole to be equal or higher to the one predicted 
by standard ACDM model, given the small value mea- 
sured by WMAP. The latter is about 5% for ACDM best 
fit quadrupoles (as can be easily estimated using a \ 2 
distribution for the quadrupole), and the difference be- 
tween the two has been the source of some confusion in 
the literature [3- 1 

The small posterior probability for ACDM model may 
still be taken as a hint that some other model may pro- 
vide a better fit to WMAP data. One could dismiss this 
result by attributing it to cosmic variance. We may sim- 
ply live in a patch of the universe where the low £ mul- 
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1 The agreement between the low qua drup ole measured by COBE 
and CDM models was discussed in llll . The better agreement 
found in that work is due to the fact that CDM models give a 
lower quadrupole than ACDM models, and that the value ob- 
served by COBE for the quadrupole was slightly higher than the 
one measured by WMAP. See also |20i'i for similar interpretations 
of the WMAP data. 
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tipoles happens to be small, for no special reason. One 
may even argue that with many independent measure- 
ments (many bins) it would be very surprising if the re- 
sults of at least one of these measurements would not 
deviate strongly from theoretical expectations. Indeed, 
WMAP shows several such peculiarities at various values 
of I. Therefore it would be premature to interpret the 
smallness of the low multipoles as a serious problem of the 
standard model, particularly in light of its great success 
in explaining and predicting major features of our uni- 
verse in great detail. Indeed this is the attitude adopted 
by Ref . Q who mention the "intriguing" smallness of the 
low I multipoles, but overall interpret their results as pro- 
viding a strong confirmation of the standard inflationary 
paradigm. However, the deficit in power observed first by 
COBE DMR and now by WMAP is an interesting result, 
and it is temp ting to look for possible explanations for 

such an effect El HI ^ ^ In this P a P er we 
discuss some candidates for the mechanism by which the 
low I multipoles may be suppressed. In other words, we 
shall investigate possible modifications of the standard 
model which can give much better probability to have 
small power in low I multipoles. 

We distinguish between two possibilities in suppressing 
the low I amplitudes. One is related to the physics of the 
inflationary phase in the early universe. Theoretically 
we have significant freedom in the design of inflationary 
potentials and, consequently, in the shape of the primor- 
dial power spectra [271 128L 12(1 Iscj. As we will show, one 
can also use the freedom in choosing initial conditions at 
the onset of inflation. This also provides a mechanism 
to affect the primordial power spectra. Both effects can 
be related through the existence of a stage when the ve- 
locity of the scalar field is not negligible, either during 
or prior to the observed 65-e-fold stage of inflation, and 
hence the perturbations cannot be described using the 
usual slow-roll parameter approximation. The drawback 
of this approach is clear; it is simply a tuning of pa- 
rameters to obtain a feature close to the present horizon 
scale. A similar problem is encountered in other attempts 

As a second option, the suppression could be due to 
late universe physics under the influence of an effective 
cosmological constant. This is probably the most inter- 
esting possibility, because it relates the suppression of 
CMB anisotropy power at horizon scales ~ H^ 1 to the 
smallness of the cosmological constant, which becomes 
dominant precisely at times ~ H^ 1 . The realization of 
such a mechanism would then link two seemingly unre- 
lated coincidence problems. 

This paper is organized as follows. In Sections [D] 
and II I II we comment on late and early universe physics 
effects as a cause for the suppression. In Section Hvl wc 
present a simple example of an inflationary model used 
to obtain a spectrum with a cutoff. In Section we fit 
two separate cutoff spectra to the WMAP observations 
in an attempt to constrain the scale of the cutoff. We 
discuss our results in Section IVll 



II. LATE UNIVERSE CASE 

Observations of late-time acceleration and low power 
on largest scales challenge our simplest theoretical ex- 
pectations. In both cases the effects become manifest on 
scales comparable to the present size of the horizon, and 
at a time equal to the present age of the universe. This 
is seen as a coincidence problem: "Why now?" . 

Indeed, CMB observations suggest that the spectrum 
of density perturbations is nearly flat on scales even a few 
times smaller than the present size of the observable uni- 
verse. Therefore, if we were born 10 9 instead of 1.4 x 10 10 
years after the big bang, we would see CMB anisotropics 
that do not display any suppression of low multipoles 
(unless we were living in an unusual part of the universe 
due to cosmic variance). Similarly, if we were living 10 9 
years after the big bang, the value of would be less 
than 1%. 

It is thus tempting to search for a physical mechanism 
relating these two problems. This may require a non- 
trivial modification of gravity at the horizon scale. For 
example, it would be appealing to have a mechanism that 
would screen not only the cosmological constant [32( , but 
also inhomogeneities on the scale 0(i/ _1 ) ~ t, on which 
the effects of the cosmological constant become mani- 
fest. Another idea is related to possible nonperturbative 
effects due to eternal inflation. For example, under cer- 
tain assumptions concerning the choice of measure in the 
theory of eternal inflation, one can come to a paradoxical 
conclusion that wc should live in a center of a spherically 
symmetric distribution of matter [3^ | . This could lead to 
suppression of the large-angle anisotropy without affect- 
ing the small-angle effects. 

However, these possibilities are very speculative. As 
a simpler alternative one may investigate whether the 
late time evolution of the metric fluctuations <& could be 
responsible for the suppression via the Integral Sachs- 
Wolfe (ISW) effect. 

The observed large-angle CMB temperature 
anisotropics are generated by scalar perturbations 
to the FRW geometry, which in the longitudinal gauge 
can be written as 

ds 2 = a 2 ( V ) ((1 + 2<f>)drj 2 - (1 - 2<S>)dx 2 ) . (1) 

The anisotropy has a contribution of ■|<I> sources at the 
last scattering time hypersurface r\ r . In addition, an inte- 
gral contribution arises, in the case where time derivative 
$' is non-vanishing. Here a{rf) denotes the scale factor, 
rj the conformal time, and a prime denotes differentia- 
tion wrt r). Consider a scale free primordial spectrum 
l^fel 2 = icr- Then, for the CMB angular power spectrum 
amplitudes we can adopt the result of |34j 

/AT\ 2 21+1 A 2 
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where the coefficient Kf is given by 



renormalizable potential 



f°° dk r 

Kf = 2£(e+l) j —[Mkrjo) 

rvo , 2 

+6 / dr]f' k j e (k(rio - 77)) 



1 



V(0) = -mV + 



(5) 



(3) 



where f' k is normalized to . The flat CDM model 
gives $' = and Kf = 1 at large scales. The flat ACDM 
model with pure cosmological constant Q\ ~ 0.7 gives 
instead iff > L amplifying the power above the Sachs- 
Wolfe plateau 0. 

However, the coefficients K] may be smaller than unity 
with particular solutions for the late time evolution of 
This can be realized if the interference between the two 
terms in the integral © is destructive, which in turn 
constrains the spectrum and evolution of the metric per- 
turbation This may require special models of the 
effective cosmological constant, or the addition of isocur- 
vature fluctuations in the quintessence field [3f| . We will 
address these possibilities in a future investigation. 



III. EARLY UNIVERSE CASE 

Here we will describe modifications of the simplest in- 
flationary models that can account for the suppression 
of the low I multipoles. As we will see, these modifi- 
cations have a common cause, which is the presence of 
a primordial or intermediate regime where some of the 
slow roll approximations arc not valid (without necessar- 
ily meaning the interruption of inflation). This results in 
a departure from the slow roll approximation when one 
evaluates the power spectrum originated during inflation. 



A. Changing the Potentials 

There have been many suggestions in the past on how 
to alter the amplitude of density perturbations produced 
during inflation, see e.g. El El El H El E3 The 
simplest idea is to fine-tune the shape of the potential. 
In general, this procedure is far from being trivial. The 
amplitude of metric perturbations at the horizon crossing 
in the slow-roll approximation is given by 



fc 3/2 $ fc = 



(4) 



where in the r.h.s. cf>(t) is a function of fc through fc ~ 
a(t)H(t), we use the units M p = (87rG)" 1/2 = 1, and V 
denotes the derivative wrt <j>. This formula is convenient 
to gain a heuristic insight of the problem. Once we start 
bending the potential, change occurs both in V(<j>) and 
in V (</>) , and the relation between <f> and the momentum 
fc of perturbations changes as well. 

To do so, one does not have to use complicated po- 
tentials. One can use, for instance, the simple generic 



By an appropriate choice of parameters we can generate 
significant features - dips or peaks - in the spectrum, as 
demonstrated in with exact solutions of the gauge 
invariant mode- by-mode equations for fluctuations. Ob- 
viously, we must tune the parameters to place the dip of 
the power spectrum around the present day cosmo- 
logical horizon. 

The situation is even simpler in the hybrid inflation 
scenario, which has more parameters by construction. 
The value of V(4>) is mainly determined by the cosmolog- 
ical constant Vq that does not change significantly during 
inflation. In this case, in order to make the amplitude of 
the potential smaller (greater) one should simply increase 
(decrease) the slope of the potential. 

Suppose, for example, that we have $ ~ 10~ 5 (as 
we should) in the theory where the potential looks like 
Vq + atfi at low 4>, corresponding to the number of c- 
folds N e < 55. Suppose also that at larger cj> the slope 
gradually increases and becomes 10 times greater on 
the scale corresponding to 60 < N e < 65. Then on 
this scale the amplitude of perturbations will become 10 
times smaller. This should lead to a primordial spec- 
trum $fc with an amplitude which is strongly suppressed 
at low fc, which will give us a desirable effect in the 
CMB anisotropy. Again, this requires some fine-tuning 
of the position of the place where the slope of the poten- 
tial changes. An analogous spectrum was found in |38j 
through a numerical mode-by-mode computation in the 
model V(<j>) ~V + A0 4 . 

Moreover, potentials of this type naturally appear in 
the simplest versions of hybrid inflation in supergravity 
(F-term inflation |3^|). If one considers N=l supergrav- 
ity with a minimal Kahler potential and a very simple 
superpotential W — S{k4>4> ~ M 2 )> one finds a potential 
of the following general form |4J| : 



V(4>)^V (l + a]n4>) + X(f> 4 + 



(6) 



Note that the potential is steep at small and at large 
</>, and relatively flat for some intermediate values of 0. 
This leads to the cutoff of the spectrum at very short 
distances, as well as on the very large scale. A similar 
behavior was advocated in 0, IE LI based on the WMAP 
data and LSS+L Q data. For certain natural values of 
parameters the infrared cutoff of the spectrum occurs on 
the scale of the horizon. This is exactly what one needs to 
explain the suppression of the low t multipoles of CMB. 



B. Kinetic regime 

It is usually assumed that the inflaton scalar field 
is moving very slowly at the beginning of inflation. 
This picture goes back to the early days of inflationary 
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paradigm, when the field was supposed to be trapped in 
the minimum of the effective potential 0,0]. Therefore 
the initial speed of the field was supposed to be zero. In 
the language of density perturbations, y ~ ^ip- ~ ^j-j 
this implied that the density perturbations produced at 
the beginning of inflation were very large, in complete 
agreement with the usual statement that the spectrum 
of density perturbations in simplest one-field inflation- 
ary models is red. 

However, in chaotic inflation the initial speed of 
the field 4> can be quite large. The field eventually slows 
down when it approaches the inflationary regime, which 
is an attractor in the phase space of all possible trajecto- 
ries {(j){t) , 4>(t)) . For a small subset of initial conditions 
the field may approach the inflationary trajectory rela- 
tively late, so that at the onset of the last 65 e-folds of 
inflation the field will have higher speed than is usually 
expected. In this case the amplitude of metric pertur- 
bations at the scale of the horizon will be smaller than 
expected, which in turn will result in a lower variance 
of CMB anisotropies at small I. This marginal inflation 
requires a tuning of the initial conditions, and not of the 
shape of the potential. 

A similar situation may occur in hybrid inflation. The 
simplest potential for two-field hybrid inflation is [16| 

V(<t>, ct) = \{a 2 - v 2 ) 2 + + \m 2 tf. (7) 

The point where <j) — <p c — and a — is a bifurcation 
point. For <f> > <f> c the squares of the effective masses of 
both fields to 2 = g 2 (f> 2 — Xv 2 + 3 Act 2 and m 2 ^ — m 2 + g 2 a 2 
are positive and the potential has a minimum at a = 0. 
For <j) < <p c the potential has a maximum at a = 0. The 
global minimum is located at (f> = and |ct| = v. Inflation 
in this model occurs while the (j) field rolls slowly from 
large values toward the bifurcation point. 

The space of initial conditions now is 
((j)(t),<j)(t);a(t),&(t)). In general evolution begins 
with large values of (f> and ct, at the boundary with 
Planck energy density. At large <p the dominant con- 
tribution to V is given by the term g 2 (j) 2 a 2 /2, so the 
Planck boundary is given by a set of four hyperbole ^3] 

<#||ct|~1. (8) 

We will assume that initially \cf>\ > \a\ and take g <C 1. 
The field <fi on this branch can take any value up to g~ x 
(for greater values of <p the field ct would have a super- 
Planckian mass gip). Suppose for defmiteness that \4>\ is 
close to its upper bound, so that the initial value of the 
field ct at the Planck boundary is of order I or somewhat 
greater. This allows for a short stage of chaotic inflation 
supported by the potential of the field ct when this field 
rolls down toward ct = 0. 

After the first stage of inflation, the field ct rapidly 
oscillates, with the frequency just slightly smaller than 
the Planck mass, and with the amplitude a(t) decreasing 



as a 3 / 2 . The oscillations induce a large contribution 
to the effective mass of the field 4>: m 2 ^ — m 2 + g 2 (<j 2 ). 
Here (ct 2 ) ~ a~ 3 stays either for the average square of 
the amplitude of the oscillations of the field, or for the 
fluctuations of this field produced during its decay to 
particles <j> and ct. 

In the beginning, g 2 (a 2 ) ^ m 2 , so the field moves with 
a much greater speed than what one could naively expect 
by looking at the effective potential of the field cj> and 
ignoring its interactions with the field ct |42l l43l| . This 
leads to a strong suppression of density perturbations 
produced at the first stages of hybrid inflation. However, 
eventually the term g 2 (a 2 } drops down and we enter the 
standard hybrid inflation regime producing perturbations 
with nearly flat spectrum. Since the term g 2 (a 2 ) drops 
down very rapidly, as a~ 3 (t), the transition between the 
flat spectrum produced at late stages of inflation and 
the strongly suppressed spectrum at the early stages of 
inflation occurs very abruptly, within a single e-fold of 
inflation. This is exactly what we want in order to explain 
the absence of the low £ multipoles. 

In terms of the slow roll parameters e = ^p- , S = — , 

one has n s — 1 — —25 — 4e. 2 Suppression of the spectrum 
at small k implies that n s — 1 , is large for k corresponding 
to the scale of the present horizon. This means that in 
our examples some of the slow-roll parameters e,<5, etc., 
at some early (or intermediate) stage of inflation were 
large. This can be attributed to the time derivatives <fi, 
</>, etc. Thus, departure from the slow roll approximation 
for calculating power spectrum is related to the regime 
where kinetic terms cannot be neglected. 

We would like to mention also, that there is a new 
possibility of modulated cosmological perturbations from 
inflation 0, where one may use a different type 
of tuning to get desirable spectra from inflation, while 
leaving the inflaton potential unchanged. 



IV. ILLUSTRATION OF THE CUTOFF 
SPECTRUM 

For illustrative purposes, we now study in a more de- 
tailed way one of the simplest possibilities we have dis- 
cussed above, namely the existence of a period of fast 
rolling of the inflaton field <fi (with strong kinetic domi- 
nation) before the last stage of inflation. The period of 
inflation is assumed to last for about 60 — 65 e-folds (the 
exact length depending on the details of reheating) such 
that the stage of fast roll can leave an imprint precisely 
at the largest scales we presently observe. As a working 



2 Note that these parameters can be further expressed in terms of 
the slow-roll parameters e, rj, A, etc., given in terms the inflaton 
potential and its derivatives. For our purposes however, it is 
convenient to keep the original parameters given in the text. 
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assumption, we start with a homogeneous and flat uni- 
verse already at the stage of fast roll, so that one may 
have to postulate the existence of a previous period of 
inflation at earlier times. Apart for this requirement, we 
leave the details of the earlier universe unspecified, and 
we rather concentrate on the evolution of the background 
and of the cosmological perturbations starting from the 
fast roll regime. 

We are mainly interested in the spectrum of primor- 
dial perturbations. We first perform an exact mode-by- 
mode numerical computation in the context of chaotic 
inflationary potential and large initial velocity for the in- 
flaton field. We will then compute the spectrum in a 
simpler idealized situation, with an instantaneous transi- 
tion between a regime of kinetic domination and a nearly 
de-Sitter stage. In this last case one can derive a simple 
analytical result which reproduces very well the exact 
spectrum obtained in the numerical evolution. 

We start with a quadratic inflaton potential, 
V ((f)) — to|</> 2 /2, and initial conditions 4>\„ = 



18.0 M p , {dcf>/dt). 



10 TOw 



giving about 61 e- 



folds of inflation. Concerning the perturbations, we work 
in momentum space and focus on the Mukhanov variable 

M 



Vk 



(9) 



where S<pk denotes the perturbations in the inflaton field, 
<E> is the metric perturbation given in Eq. Q), prime de- 
notes derivative with respect to conformal time rj, and 
h = a! I a , a being the scale factor of the universe. The 
evolution of is given by (see e.g. 01 f° r details) 

, 9 z"~\ ad>' 

k 2 v = , z=-f- (10) 

z h 



v" + 



Due to the choice of initial conditions, the system is 
initially in a kinetic dominated regime, with h oc <fi' . As 
a consequence, we have initially 



a 



2h rj 



hi 



{l + 2hoT)Y 



(11) 



with ho denoting the value of h at the time rj — 0, where 
we have normalized a = 1 . 

As long as the approximation l|ll|) is valid, Eq. (|10fl is 
solved by 



Vk (v) 



8 h 



k 
2h~o 



(12) 



(2") 

where Hq is a specific Hankel's function and where we 
have chosen (and appropriately normalized) the solution 
of Eq. HlOfl which reduces to exp (— i krf) j\/2 k at very 
short wavelengths, k ^S> 1/M • 3 



The choice {HJ can be motivated by taking an adiabatic vacuum 



In the numerical integration we have solved exactly the 
evolution equations for the background quantities 4> & n d 
a , together with Eq. (fTTTfl for the perturbation mode Vk ■ 
We have used Eq. I|12l) only to set the value of Vk and its 
derivative at the initial time rj = . We have plotted in 
the last panel of fig. [21 the final spectrum for the metric 
perturbation . 

Modes with very large k are well inside the horizon 
during the fast roll of 4> j and hence they are quite insen- 
sitive to the evolution of the background at this stage. 
Their horizon crossing occurs during inflation, and the 
usual computation (giving nearly scale invariant spec- 
trum) applies. At wavelengths comparable with the size 
of the horizon H~ x during the onset of inflation we find 
instead a strong departure from scale invariance. The 
amplitude of the spectrum is oscillating for k >, , fol- 
lowed by a sharp decrease at k ~ H* . By an appropriate 
choice of the number of e-folds of inflation, it is possible 
to relate such wavelengths to the size of the present hori- 
zon, with a consequent decrease of the observed CMB 
anisotropics at large scales. 

This result is well reproduced in a highly simplified 
situation, with an instantaneous transition between the 
regimes of kinetic domination and nearly de-Sitter expan- 
sion. The scale factor a in the idealized model evolves as 



a ~ ^l + 2Hi 1 , 7] < 

a ~ - — l -— , 7] > 
1 — H 7] 



(13) 



where we have now set rj = , and a = 1 at the transition, 
while H denotes the (physical) Hubble parameter during 
inflation. 

The evolution equation for the Mukhanov variable Vk 
is again given by Eq. I|10[) . with z"/z ~ a" /a in both 
regimes. During kinetic domination, the evolution of Vk 
is again given by Eq. I|12|) with the quantity ho replaced 
by H . During the de-Sitter stage one finds instead 



Ce -i(k v -k/H) I j _ 



kr}- k/H 



kr] — k/H 



(14) 



and the coefficients of the two modes can be obtained by 
requiring continuity of v and v' at the transition, 



G 



-ik/H 



a/32 H/tt 



H 



(2) 



2H 



during a previous inflationary stage, and assuming an adiabatic 
transition into the kinetic dominated regime (adiabaticity may 
be preserved due to a high effective mass for the field during 
these previous times, as for example in the hybrid inflationary 
case described in the previous secti on) . An analogous initial 
state was considered for example in l48l . where the production 
of gravitational waves was studied in a model with a radiation 
dominated stage followed by inflationary expansion. 



6 



1.2 



K 08 
I 0.6 

a 



0.4 



0.2 




10 



k/H 



FIG. 1: Power spectrum for Q from the analytical compu- 
tation, Eqs. 11151 and 11611 . This is to be compared with the 
spectrum of $ obtained by the numerical evolution, and re- 
ported in the last panel of fig. 



D 



ik/H 



^32 H/n 



H, 



(2) 



k 
2H 



k 1 \2H 



(15) 



We consider the spectrum of Q = v/a which becomes 
constant at late time, as indicated by the leading order 
contribution to equation l|14fl for i] — ► 1/H 



P Q =*\Q\*^f- k \C-D\* 



(16) 



At short wavelengths (fc 3> H) one recovers the standard 
result of a (nearly) scale invariant spectrum 



\C\ 

Pq 




\D\ « \C\ 



(17) 



At super-horizon scales, the two modes Qk and are 
related by a k independent rescaling so that the spectrum 
given by Eqs. (|15fl and (jlfij) directly translates into the 
spectrum of $ , up to an overall normalization factor. 

The spectrum obtained by this simple analytical cal- 
culation is shown in figure ^ We see that it reproduces 
very well the spectrum of $ obtained with the exact nu- 
merical evolution, which is reported in the last panel of 
fig- 121 



V. CMB ANISOTROPIES WITH A CUTOFF 

To determine how strongly the data favour a primor- 
dial power spectrum with a feature at large scales we 
have fit the WMAP data with two template cutoffs. The 
first is derived from the simple model discussed above 



and includes the detailed oscillatory effect at the tran- 
sition between the two regimes in the initial spectrum. 
The second is based on an exponential cutoff of the form 



fc 3 |$(fc)| 2 = A s (l -exp 



-(fe/fec) c 



)k n " 



(18) 



where we have chosen a to fit the rise of the model with 
a = 3.35. A s is the overall amplitude of the scalar pertur- 
bations. This parametrization does not include features 
in the power spectrum which are dependent on the spe- 
cific model considered, but which are too small to change 
qualitative features in the CMB power spectrum. 

Due to the nature of the TE cross-correlation signal re- 
ported by WMAP , we chose to fit our models to the TT 
data only. We believe that it is still premature to draw 
too strong conclusions by making use of the E-type po- 
larization signal, which has only been measured through 
its correlation with the much larger temperature signal 
by WMAP. Once the WMAP E-type observations have 
reached a sufficiently high signal-to-noise level to unam- 
biguously measure a polarization signal free of systemat- 
ics, its use in constraining models will be much safer. 4 

We use a modified version of the publicly available 
CMBFAST |H3 code to compute CMB spectra with vari- 
ous k c in the two cutoffs. We generate grids of models by 
varying three parameters which affect the shape of the 
spectrum at the lowest multipoles, namely the spectral 
index of the primordial power spectrum n Sl the energy 
density of cosmological constant component in units of 
the critical density and the present day wavenum- 
ber of the cutoff in the spectrum k c in units of Mpc -1 . 
The grid values are regular in the parameters and have 
a resolution of 32, 20, and 32 grid points with ranges 
[0.83, 1.04], [0.6, 0.85], and [2.0 x 10~ 3 , 1.0 x 10~ 5 ] respec- 
tively. 

At each point in the grid we use subroutines derived 
from those made available by the WMAP team to eval- 
uate the log likelihood with respect to the WMAP data 
[5l| . Other parameters that are not expected to affect 
the shape of the spectrum at I < 30 are held fixed at 
the following values; f^/i 2 = 0.022, tt m h 2 — 0.135 and 
t c = 0.17 (3- In addition we consider only flat models 
with £! t ot = 1- The present value of the Hubble con- 
stant for each model is then fixed for each choice of f^A 
in the grid. As we are not interested in the overall ampli- 
tude of the primordial perturbations we marginalize over 



4 An analysis with inclusion of TE data was first performed in l2 r i , 
and, successivel y to li r l and to the firs t version of the present 
manuscript, in |4!.I . The results of |4'1 are similar to ours, but 
their conclusions concerning the problem of low I multiples are 
slightly more optimistic due to the strong emphasis put on the 
TE data. We caution that taking into account both cosmic vari- 
ance and the correlation with the much larger TT signal would 
result in a signifi can t increase of the error bars on the TE data 
(see fig. 2 of ref . l25l , where this is explicitly shown for a partic- 
ular model), possibly affecting claims which strongly depend on 
their use. 
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a range of values in the amplitude of the power spectrum. 
We only consider scalar perturbations in this work. 

Our choice of parameters is motivated by our interest 
in motivating the requirement for a cutoff in the spec- 
trum and not by an attempt at a precise determination 
of cosmological parameters, an exercise which should be 
the focus of much more exhaustive investigations 0, |25| . 
Residual correlations with parameters that have been 
fixed such as r c will be subsumed by our tilt and am- 
plitude parameters which should be viewed as mildly ef- 
fective quantities. In addition our use of a prior on Q\ 
(see below) is motivated by a simple mimicking of Large 
Scale Structure (LSS) and Type-la supernovae (SNla) 
constraints external to the WMAP results. 5 

In Fig. |2 we show the resulting best-fit models we ob- 
tain from the grid for both cutoff models considered. The 
top panel shows a standard (no cutoff) spectrum model 
with the same parameters as those obtained in the best- 
fit for the two models (n s = 0.9587 and Q A = 0.732). We 
have displayed the noise contribution to the uncertainties 
as errors on the binned WMAP results. The cosmic vari- 
ance contribution to the errors is shown as upper and 
lower contours around the model. We have displayed 
the la confidence limits for a lognormal distribution as 
this approximates more clearly the significance of the low 
quadrupole and octupole signal. The middle two panels 
show the best-fit for the model cutoffs. We can see that 
both models fit the lower multipoles better then the stan- 
dard spectrum by reducing power on the largest scales. 
The bottom panel shows the two template spectra used 
in the fits shifted to an arbitrary cutoff scale of k c — 1 
for a n s = 1 model. 

The bottom panel shows the two template spectra used 
in the fits shifted to an arbitrary cutoff scale of k c = 1 for 
a,n s — I model. It is important to note that even a step- 
like cutoff in the initial spectrum will not produce a step- 
like feature in the CMB spectrum due to the convolution 
of the initial spectrum in Eq. [21 Thus it is difficult to 
reproduce the sharp drop-off that occurs in the angular 
power spectrum. 

Although we fit for all TT multipoles in the range 
2 < £ < 900, most of these simply provide an accu- 
rate calibration of the overall amplitude of the model 
as we are keeping the matter content fixed. We therefore 
concentrate on the lowest multipoles in considering the 
significance of our best-fit models. For £ < 10 we find 
that the probability of exceeding the observed \ 2 lS 5.6% 
in the case of the standard model, 21.7% when using the 
template model cutoff, and 33.4% with the exponential 



6000 



4000 



5 The authors of l4!l erroneously compare our analysis to their 
"WMAP TT only" results. More appropriately, our simple 
choice of parameters and priors is closer to a "WMAPext TT + 
2dF + SNla" combination. The latter requires a very different 
parameter grid compared to the "WMAP TT only" data, whose 
best fit model lies out in the tail of the LSS+SNla motivated 
priors 
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k/(aH). 



FIG. 2: CMB anisotropies with cutoff primordial spectra. 
The top panel shows a standard HaCDM model without cut- 
off, and with n s = 0.9587 and Qa = 0.7316, corresponding to 
the best-fit values of our model grids. The upper and lower 
contour show the la confidence level for a lognormal distri- 
bution appropriate for the sample variance of the observa- 
tions. The error bars show the experimental noise contribu- 
tions. The result in the second panel is obtained with the 
same model, but with an infrared cutoff on the primordial 
power spectrum, as obtained from an exact evolution of the 
model described in section lTVl The result in the third panel is 
instead obtained with the power spectrum given in Eq. 1181 . 
The two cutoffs chosen correspond to the best-fit cutoff scales 
k c = 4.9xl0" 4 Mpc _1 and k c = 5.3 x 10" 4 Mpc _1 respectively. 
The bottom panel shows the two primordial power spectra 
corresponding to these cutoffs and for the case n 3 = 1. 



cutoff of Eq.UHl 

In Fig. El we marginalize over the range in fl\ with a 
Gaussian prior fl\ — 0.70 ± 0.03 and plot the resulting 
X 2 values as a function of n s and cutoff scale k c . The 
contours shown are for A\ 2 values giving one, two, and 
three a contours for two parameter Gaussian distribu- 
tions. We find that the distribution for the cutoff scale 
has a tail as k c — ► close to the 2a level. This is evident 
from the fact that the models asymptotically approach 
the scale invariant case which has a low but not vanish- 
ing probability with respect to the observations [3]. The 
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FIG. 3: Marginalized distribution in the n a -k c plane for our 
grid of models. The left panel is for the exponential cutoff 
while the right panel is for the toy model discussed in the text. 
The solid contours show drops in \ 2 corresponding to 68%, 
95%, and 99.7% confidence for a two-dimensional Gaussian 
distribution. 
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FIG. 4: The one dimensional marginalized distribution for 
the cutoff scale k c . The solid (red) curve is for the exponen- 
tial cutoff while the dashed (blue) curve is for the spectrum 
obtained in the inflation model discussed in the text. Notice 
that the sharp drop-off in the likelihood at high k c is driven 
mainly by our prior on Oa- 



data however prefer a cutoff scale within the present ob- 
servable universe at the la level. 

Finally we marginalize over n s with a flat prior to ob- 
tain the marginalized probability distribution in k c shown 
in Fig. |Ij As the distributions do not vanish as k c — > we 
integrate the functions in the range < k c < oo to obtain 
68% and 95% confidence limits. For the inflation model 



we find k c = 5.3±£{j x lO^Mpc" 1 with a 95% upper limit 
of k c < 7.02 x 10 _4 Mpc _1 . Similarly for the exponen- 
tial cutoff model we find k c = 4.9^^ x 10 _4 Mpc~ 1 with 
a 95% upper limit of k c < 7.4 x lO^Mpc" 1 . It is im- 
portant to note that the confidence limits depend on the 
integration measure adopted. In this case we have taken 
a measure linear in k c ; however an alternative would be 
to take a measure in In k c which may be more suitable 
for the distribution of the cutoff scales given by inflation 
theories. This second option, which gives higher weights 
to lower wavenumbers, will crucially depend on the value 
chosen for the largest scales, In A:™ 111 > — oo. In the con- 
text of inflationary models it is natural to limit to those 
scales which left the horizon when the energy density of 
the universe was close to Planckian. 

When faced with such ambiguities it is useful to refer 
to the difference in likelihoods between different values 
of k c . This allows an approximate estimate of the sig- 
nificance of having a cutoff without depending on any 
measure used in deriving confidence limits. As an exam- 
ple the difference in \ 2 between the peak in the distri- 
butions and at k c = is found to be A% 2 = 5.3 for the 
exponential cutoff and A% 2 = 4.2 for the inflation model. 
This indicates that the case with no cutoff is just over 2a 
below our best-fit cutoff model. 



VI. DISCUSSION 

Both COBE DMR and WMAP observations indicate 
a low value for the power on largest scales. Interest- 
ingly, these scales are of the same order as the scale cor- 
responding to the time of dominance of the cosmological 
constant. If one ignores the possible interpretation of 
the suppression of low multipoles due to cosmic variance 
(which is perhaps the simplest and the least dramatic 
interpretation), one faces a new coincidence problem in 
addition to the usual one required to explain the late- 
time acceleration of the universe (i.e. Q\ ~ 20 m today). 
In this paper we briefly discussed the possibility that 
the late-time evolution of the metric perturbations 
may reduce the power in the low multipoles through inte- 
grated Sachs- Wolfe effect. This may require a combined 
investigation of adiabatic and isocurvature perturbations 
in certain models of dark energy. 

The main emphasis of our paper was related to mod- 
ifications of the simplest models of inflation due to the 
possible existence of a kinetic stage, when the velocity 
of the scalar field was not negligible. The breaking of 
the slow roll approximation which could occur at this 
stage may cause a significant suppression of the large 
scale density perturbations. In order to affect the low 
£ multipoles, this stage should occur in the beginning of 
the last 65 e-fold period of inflation. Clearly, this requires 
either fine tuning of parameters in the inflaton potential, 
or fine-tuning of initial conditions which are necessary to 
obtain a desirable modification of the spectrum close to 
the present horizon scale. 
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The final part of the paper has the aim to quantify 
with a simple example the level of suppression required 
by the data. A stage of kinetic domination before the 
observable stage of inflation can result in a sharp cutoff 
in the spectrum of primordial perturbations. 

Fits to the observations reveal that the data favour a 
cutoff in the spectrum although the significance of the 
"detection" is less than 3er. This is a reflection of the 
fact that the k c = case (i.e. standard ACDM model) 
is not a terrible fit to the data. The price of obtaining 
a better fit is the addition of an extra parameter in the 
theory which should probably be done under stronger 
observational evidence than we currently have. 

On the other hand, we should note that when we at- 
tempt to evaluate all possible theories that could produce 
the required suppression of the spectrum, we consider 
the space of all possible cutoff parameters equally dis- 
tributed in momentum space. Meanwhile if we would 
assume equal prior distribution in space of parameters of 
inflationary models, this would translate into an approx- 
imately equal distribution in terms of In k rather than 
k. This remark also applies to other attempts to explain 
the suppression of low multipoles 0, IH, 0] and 
constitutes a choice of prior on k c a problem encountered 
when constraining many of the cosmological parameters 
(e.g. a choice of uniform prior on Hq as opposed to in 
J7a)- However, as we have discussed, a possibly less am- 
biguous quantity is the likelihood ratio between different 
values of the cutoff since this does not depend on any 
underlying theoretical distribution. We have quantified 
this in the previous section. 

We should also bear in mind that in any case, observa- 
tions tell us the values of the parameters, and we should 
accept them even if they seem unusually fine-tuned. For 
example, many coupling constants in the standard the- 
ory of electroweak interactions are of the order lCP 1 . 
Meanwhile the coupling constant of the electron and the 
Higgs field, which is responsible for the electron mass, 
is 2 x 10~ 6 . The smallness of this parameter is notable, 
but usually we do not assume that all coupling constants 
are random variables and do not calculate the probabil- 
ity to have such a small coupling constant on the basis of 
its comparison to other coupling constants in the theory. 
We still do not know why this coupling constant is so 



incredibly small, but it does not make the standard the- 
ory of electroweak interactions unacceptable. Moreover, 
according to the standard model, the neutrino must be 
massless, meanwhile the recent data suggest that this is 
not the case. This means that we need to extend the 
standard model by making it even more complicated. 
However, this does not cast any doubts concerning the 
basic idea of gauge invariance and spontaneous symmetry 
breaking, which made the electroweak theory internally 
consistent. 

High energy physicists never claim that they fine-tune 
the parameters of the theory, they simply fit the data. 
Until very recently cosmologists did not need to do the 
same and had the luxury to debate which values of the 
parameters seemed more natural. We are now entering 
the age of precision cosmology, and we need to adjust our 
attitude accordingly. 

During the last 20 years, the inflationary scenario 
emerged as the simplest and possibly even unique way 
of constructing an internally consistent cosmological the- 
ory. The main result of our paper is that, if needed, 
we can fit the parameters of the inflationary theory and 
make it consistent with the small amplitude of the large- 
angle anisotropy of the CMB. It is incomparably easier 
to do so than to explain all features of the spectrum of 
the CMB, as well as the observed homogeneity, isotropy 
and flatness of the universe without using inflation. 
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